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EXECUTIVE  SUMMARY 


A  critical  issue  identified  by  the  aviation  industry  is  the  need  to  examine  the  effects  of  repairs  on 
the  structural  integrity  of  aircraft.  The  incorporation  of  damage  tolerance  methodologies  in  the 
maintenance  and  repair  practices  of  aging  aircraft  is  required  in  order  to  insure  their  continued 
airworthiness  and  operational  safety.  In  December  1978,  the  Federal  Aviation  Administration 
(FAA)  amended  their  Fatigue  Evaluation  Requirements  to  include  a  damage  tolerance  philosophy. 
However,  the  majority  of  the  current  aircraft  repairs  are  still  designed  using  static  strength 
approaches.  The  resources  needed  for  damage  tolerance  designs  of  repairs  are  lacking, 
particularly  for  small  operators  and  independent  repair  facilities.  Consequently,  inadequate  repairs 
are  being  designed  that  do  not  meet  damage  tolerance  requirements.  In  an  effort  to  address  this 
need,  a  task  was  undertaken  to  develop  a  new  user-fiiendly  software  tool.  Repair  Assessment 
Procedure  and  Integrated  Design  (RAPID),  capable  of  static  strength  and  damage  tolerance 
analyses  of  simple  fuselage  skin  repairs. 

A  simplified  engineering  approach  to  static  and  damage  tolerance  analyses  of  riveted  fuselage  skin 
repairs  has  been  incorporated  in  RAPID.  In  this  study,  the  damage  tolerance  analysis 
methodology  in  RAPID  was  evaluated  in  terms  of  the  fastener  loads,  stress-intensity  factor 
solutions,  crack  growth,  residual  strength,  and  inspection  schedule  calculations.  Three  example 
problems,  each  representing  a  typical  fuselage  skin  repair  configuration,  were  analyzed.  The 
analysis  results  obtmned  fi'om  RAPID  were  compared  with  results  generated  using  a 
Representative  Original  Equipment  Manufacturer  (ROEM)  method  and  a  special  purpose  finite 
element  program  for  fiacture  mechanics  analysis  and  crack  growth  simulation  in  layered  two- 
dimensional  structures.  In  general,  results  generated  using  RAPID  were  in  good  agreement  with 
results  generated  using  the  ROEM  method  and  the  finite  element  code. 
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INTRODUCTION 


A  critical  issue  identified  by  the  aviation  industry  is  the  need  to  examine  the  effects  of  repairs  on 
the  structural  integrity  of  aircraft.  The  incorporation  of  damage  tolerance  methodologies  in  the 
maintenance  and  repair  practices  of  aging  aircraft  is  required  in  order  to  insure  their  continued 
airworthiness  and  operational  safety.  In  December  1978,  the  Federal  Aviation  Administration 
(FAA)  amended  their  Fatigue  Evaluation  Requirements  to  include  a  damage  tolerance  philosophy. 
However,  the  majority  of  the  current  aircraft  repairs  are  still  designed  using  static  strength 
approaches.  The  resources  needed  for  damage  tolerance  designs  of  repairs  are  lacking, 
particularly  for  small  operators  and  independent  repair  facilities.  Consequently,  inadequate  repairs 
are  being  designed  that  do  not  meet  damage  tolerance  requirements. 

In  an  effort  to  address  industries  need  for  damage  tolerance  repair  analysis  capabilities,  a  task  was 
undertaken  to  develop  a  new  user-friendly  software  tool.  Repair  Assessment  Procedure  and 
Integrated  Design  (RAPID).  The  software  tool  conducts  static  strength  and  damage  tolerance 
analyses  of  simple  fuselage  skin  repairs  and  can  be  used  to  assess  repair  design  and  integrity, 
provide  advice  on  repair  design  improvement,  and  provide  data  to  establish  an  inspection  program 
for  the  repair.  The  system  is  designed  so  that  small  independent  repair  stations  and  operators  of 
small  commuter  airplanes  will  be  able  to  use  the  tool. 

RAPID  is  a  simple,  PC-based  repair  tool  with  a  modularized  open  system  architecture  providing 
easy  implementation  of  tool  upgrades  and  new  features.  The  current  version  of  RAPID  consists 
of  five  modules:  Graphics  User  Interface  (GUI),  Advisory  System  (AS),  Static  Analysis  (SA), 
Damage  Tolerance  Analysis  (DTA),  and  Database  modules.  The  GUI  module  is  Windows"™- 
based  providing  user-friendly,  point-and-click,  multitasking  capabilities.  The  AS  module  provides 
generic  engineering  guidelines  for  specifying  fastener  arrangements,  doubler  materials,  and 
doubler  thicknesses.  The  SA  module  performs  static  strength  analysis  of  repair  configurations 
based  on  the  fastener  shear  allowables,  skin  hole  bearing  allowables,  and  doubler  tensile 
allowables.  The  DTA  module  calculates  the  critical  fastener  loads  and  applies  a  simplified  crack 
growth  algorithm  to  determine  crack  growth  and  residual  strength  and  provides  an  inspection 
schedule.  The  Database  module  consists  of  material  database  conteuning  elastic,  strength,  crack 
growth,  fi’acture  properties,  and  a  fastener  database  containing  fastener  shear  allowables  and  skin 
hole  bearing  allowables. 

In  this  study,  the  damage  tolerance  analysis  methodology  in  RAPID  was  evaluated  in  terms  of  the 
fastener  loads,  stress-intensity  factors,  crack  growth,  residual  strength,  and  inspection  schedule. 
Three  example  problems,  each  representing  a  typical  fuselage  skin  repair  configuration,  were 
analyzed.  The  analysis  results  obtained  from  RAPID  were  compared  vrith  results  generated  using 
a  Representative  Original  Equipment  Manufacturer  (ROEM)  approach  and  the  Fracture  Analysis 
Code  for  2-Dimensional  Layered  (FRANC2DL)  structures,  a  special  purpose  finite  element 
program  for  fracture  mechanics  analysis  and  crack  growth  simulation  in  layered  two-dimensional 
structures. 

In  the  following  sections,  the  damage  tolerance  methodology  is  highlighted  and  the  models  used 
in  the  FRANC2DL  and  the  ROEM  analysis  are  presented.  The  results  are  then  compared  in  terms 
of  the  fastener  loads,  stress-intensity  factor  solutions,  crack  growth  characteristics,  residual 
strength,  and  inspection  schedules. 
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APPROACH 


RAPID  DAMAGE  TOLERANCE  METHODOLOGY. 


The  damage  tolerance  analysis  in  RAPID  is  based  on  the  approach  outlined  by  Swift  [1].  The 
procedure  involves  several  steps  including  the  calculation  of  fastener  loads,  assumption  of  the 
initial  flaw  geometry  and  growth,  and  calculations  of  stress-intensity  factor,  crack  growth, 
residual  strength,  and  inspection  schedule.  Each  of  these  components  is  discussed  in  the 
foUowdng  paragraphs. 

FASTENER  LOADS.  A  single-strip  model  is  used  to  determine  the  stresses  in  the  vicinity  of  the 
critical  fastener  hole.  In  the  analysis,  the  skin  bypass  stress  skin  bearing  stress 

^bearin^  at  the  critical  fastener  hole  are  calculated  due  to  the  applied  tensile  gross  stress  (Jg^oss  • 
The  skin  and  doubler  are  idealized  as  axial  rods  having  stififhesses  equal  to  AE/L  where  A  is  the 
cross  section  area,  E  is  the  modulus,  and  L  is  axial  length.  The  fastener  stiffness  is  expressed 
empirically  as  [1]: 


^fast  “■ 


E'D 


A  +  B\ 


(1) 


where  A  and  B  are  curve  fitting  parameters,  E'  is  the  effective  modulus  of  the  skin  and  doublers, 
and  are  the  skin  and  doubler  thicknesses,  and  D  is  the  fastener  diameter.  For  aluminum 
fasteners,  A  =  5.0  and  B  =  O.S.  Using  a  standard  matrix  method  of  structural  analysis,  the  skin 
internal  loads,  crijypas^  and  (Jbearing  >  can  be  determined  and  expressed  in  terms  of  • 

^bypass  ~  gross  >  ^bearing  ~  gross  (2) 

where  ^is  the  ratio  of  bypass  to  gross  stress,  and  ^is  the  ratio  of  bearing  to  gross  stress. 

INITIAL  FLAW  GEOMETRY  AND  GROWTH.  An  initial  flaw  geometry  is  assumed.  The  two 
initial  crack  cases  available  in  RAPID  are  shown  in  figure  1.  For  the  first  case,  the  initial  flaw 
consists  of  a  primary  crack,  Op  =  0.05",  and  a  secondary  crack,  Og  =  0.005",  at  the  edge  of  the 
center  fastener  hole.  For  the  second  case,  the  initial  flaw  consists  of  a  single  crack,  Up  =  0.05",  at 
the  edge  of  the  center  fastener  hole.  The  subsequent  crack  growth  scenarios  are  also  illustrated  in 
figure  1.  For  the  first  case,  when  the  primary  crack  grows  into  the  first  adjacent  hole,  a  crack  of 
length  0.005"  is  then  assumed  to  initiate  on  the  opposite  side.  For  the  second  case,  when  the 
primary  crack  grows  into  the  adjacent  hole,  cracks  of  length  0.005"  are  assumed  on  the  opposite 
side  of  each  of  the  two  holes.  Whenever  a  crack  enters  a  hole,  a  crack  of  length  0.005"  is 
assumed  to  initiate  on  the  opposite  side  of  the  hole. 
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Initial  Flaw 


a,  =  0.005"  ap  =  0,05" 

c.„.o  oWo  o 

Subsequent  Growth 

Crack  Grows  '  a  =  0.005" 

Into  Hole 


-a  =0.05" 


Case  2 


O  O  0^0  o  Initial  Flaw 

CZ^  Subsequent  Growth 
a  =  0.005"  Crack  Grows  a  =  0.005" 

Into  Hole 


FIGURE  1.  THREE  INITIAL  FLAW  GEOMETRIES  AND  SUBSEQUENT  GROWTH 

STRESS-INTENSITY  FACTOR.  There  are  five  baseline  stress-intensity  factor  (SIF)  solutions  in 
RAPID  as  shown  m  figure  2:  (1)  a  single  crack  emanating  fi-om  a  hole  in  an  infinite  plate 
subjected  to  far-field  tension;  (2)  a  single  crack  emanating  fi-om  a  hole  in  an  infinite  plate 
subjected  to  a  pair  of  pin  loads;  (3)  two  cracks  of  unequal  length  emanating  fi-om  a  hole  in  an 
infinite  plate  subjected  to  far-field  tension;  (4)  two  cracks  of  unequal  length  emanating  fi-om  a 
hole  in  an  infinite  plate  subjected  to  a  pair  of  pin  loads;  and  (5)  a  crack  in  an  infinite  plate 
approaching  an  open  hole  subjected  to  far-field  tension.  Methods  of  superposition,  compounding, 
and  similarity  are  used  to  determine  the  SIF  needed  during  the  simulation  of  crack  growth.  The 
SIF  is  expressed  in  terms  of  the  applied  stress,  S,  crack  length,  a,  and  a  boundary  correction 
factor,  p,  as  K  = 


/  \ 

/ 

;  ; 

O 

j 

"r 

V 

r 

2. 

4. 

( r 

FIGURE  2.  FIVE  BASELINE  STRESS-INTENSITY  FACTOR  SOLUTIONS  IN  RAPE) 
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CRACK  GROWTH  MODEL.  A  simplified  crack  growth  model  based  on  Walker’s  equation  is 
used; 

where  R  is  the  stress  ratio,  K^ax  is  the  maximum  stress-intensity  factor  (due  to  CTbypass  >  ^bearing  > 
and  and  C,  p,  and  q  are  material  coefficients.  Inverting  and  integrating  equation  3  and 

separating  geometry,  load,  and  material  terms  3nelds  an  expression  for  the  number  of  cycles. 


The  first  terms  in  front  of  the  integral  are  the  load  and  material  terms  where  S  is  the  maximum 
applied  stress.  In  this  study,  a  constant  amplitude  stress  is  assumed  where  R  =  0mdS=  (Tgross- 
The  integral  is  the  geometry  term  where  a,-  is  the  initial  crack  length,  ay  is  the  final  crack  length, 
and  y5is  the  boundary  correction  factor.  Gauss  numerical  integration  is  used  to  solve  equation  4. 

RESIDUAL  STRENGTH.  The  residual  strength  is  then  calculated  in  terms  of  the  skin  plane 
stress  fracture  toughness  the  crack  length  a,  and  the  boundary  correction  factor,  /3. 

(5) 

The  inspection  threshold  is  calculated  in 


c 

i^res  ~  I — 

INSPECTION  THRESHOLD  AMD  INTERVAL. 
RAPID  from 


where  iVc„/is  the  number  of  flights  for  the  crack  to  propagate  from  an  initial  flaw  size  of  0.05"  to 
the  critical  size,  cicrit’  under  the  limit  load  condition.  The  limit  load  condition  was  derived  from 
the  damage  tolerance  evaluation  requirements  specified  in  FAR  §25.571  and  is  assumed  to  be  the 
following; 


Ll{P  +  0.5)r 


where  P  is  the  pressure  differential  of  the  fuselage,  r  is  the  fuselage  radius,  and  t  is  the  fuselage 
skin  thickness.  The  critical  crack  length  is  given  in  terms  of  the  limit  load  condition  and  plane 
stress  fracture  toughness; 


The  inspection  interval  is  calculated  in  RAPID  by  the  following  equation; 

yj  _  ^crit  ~  ^det 

ly  ivit  —  : 
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=  157,500  Ib/in  calculated  from  equation  1.  The  fastener  arrangement  consists  of  five  rivet 
rows  per  repair  side  with  a  pitch  distance  of  1"  and  an  edge  distance  of  0.5".  An  initial  flaw  was 
assumed  at  the  edge  of  the  middle  hole  in  the  first  fastener  row  consisting  of  a  primary  crack  of 
length  Cp  -  0.05"  and  a  secondary  crack  of  length  =  0.005".  A  remote  tensile  stress  of  ~ 
17  ksi  was  applied  to  the  skin. 

Two  idealizations  shown  in  figure  4  were  used  to  represent  the  repair  configuration  shown  in 
figure  3.  First,  a  single-strip  model  was  used  to  determine  the  bearing  stress,  <ybearing  » 
bypass  stress,  bypass  >  the  critical  fastener  hole.  Then  a  plate  (skin)  model  containing  a 
periodic  array  of  holes  was  used  to  determine  the  SIF  due  to  the  applied,  bearing,  and  bypass 
stresses.  For  the  loading  on  the  single-strip  model,  the  results  from  RAPID  were  used  where 
Ogross  =17  ksi,  cTbearing  ~  27.0667  ksi,  and  cTbypass  =11-92  ksi.  These  stresses  were  then  applied 
to  the  plate  model  and  the  SIF  was  calculated  at  various  crack  lengths  up  to  the  second 
breakthrough  condition.  Breakthrough  is  defined  as  the  condition  where  the  crack  grows  into  the 
next  hole.  Symmetric  loading  conditions  were  applied  to  the  plate  model  using  the  method  of 
superposition. 


Single-Strip  Model  Plate  Model 

RAPID  FRANC2DL 


FIGURE  4.  MODELS  USED  TO  REPRESENT  REPAIR  CONFIGURATION 

Three  cases  with  various  loading  and  fastener  arrangements  were  analyzed  for  the  plate  model  as 
shown  in  figure  5.  For  Case  1,  a  plate  containing  nine  fastener  holes  was  subjected  to  a  remote 
stress  of  l/2(crg^oj^  +  cxbypas^  a  pair  of  bearing  loads  y^iPbectringt  ^1  fastener  holes.  For 
Case  2,  a  plate  containing  nine  fastener  holes  was  subjected  to  two  loading  conditions:  (1)  from 
the  initial  flaw  to  the  first  breakthrough,  a  remote  tensile  stress  of  l/2(c^05,^  +  %ypas^ 
applied  in  the  far-field  and  a  pair  of  bearing  loads  ^^{(Jbearin^  was  applied  to  the  middle  fastener 
hole 
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Case  2 


Initial  Flaw  To  First  Breakthrough  First  To  Second  Breakthrough 


FIGURE  5.  APPLIED  LOADING  AND  FASTENER  ARRANGEMENT  FOR  THE 

FRANC2DL  MODELS 


containing  the  initial  flaw,  and  (2)  from  the  first  breakthrough  to  the  second  breakthrough,  a 
remote  tensile  stress  of  l/2(o^^o^^  +  cri^yp^^)  was  applied  in  the  far-field  and  a  pair  of  bearing 
loads  l/2(<Ji,garing)  ^as  applied  to  the  middle  three  fastener  holes  which  were  connected  by 
cracks.  For  Case  3  two  different  fastener  and  loading  arrangements  were  considered:  (1)  from 
the  initial  flaw  to  the  first  breakthrough,  a  plate  containing  three  fastener  holes  was  subjected  to  a 
remote  tensile  stress  of  l/2(agross  +  <^bypass)  ^  of  bearing  loads  '^l'2.{aijearin^ 
middle  fastener  hole  containing  the  initial  flaw,  and  (2)  from  the  first  breakthrough  to  the  second 
breakthrough,  a  plate  containing  five  fastener  holes  was  subjected  to  remote  tensile  stress  of 
l/2(cj^o_y^  +  CTijypas^  and  a  pair  of  bearing  loads  y^(,<^bearin^  middle  three  fastener  holes 
which  were  connected  by  cracks.  Of  these  three  cases.  Case  3  most  closely  resembles  the  loading 
and  fastener  arrangement  conditions  used  in  the  RAPID  damage  tolerance  analysis  described 
earlier. 

Figure  6  shows  the  two  finite  element  meshes  which  were  used  to  analyze  the  three  cases  studied. 
The  first  mesh  was  used  for  crack  lengths  between  the  initial  flaw  up  to  the  first  breakthrough 
condition.  This  mesh  consisted  of  3008  eight-noded  elements.  The  second  mesh  was  used  for 
crack  lengths  from  the  first  breakthrough  up  to  the  second  breakthrough  condition.  This  mesh 
consisted  of  3172  eight-noded  elements.  Both  meshes  changed  automatically  during  the  analysis 
due  to  the  adaptive  meshing  used  in  FRANC2DL  when  cracks  are  introduced. 


Initial  Flaw  To  First  Breakthrough 


Detailed  View: 

3008  Eight-Noded  Elements 


First  To  Second  Breakthrough 


Detailed  View: 

3172  Eight-Noded  Elements 


FIGURE  6.  FINITE  ELEMENT  MESHES  USED  IN  THE  FRANC2DL  ANALYSIS 
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CRACK  GROWTH  AND  RESIDUAL  STRENGTH.  Eight  steps  were  used  to  calculate  the 
crack  growth  and  residual  strength: 

1 .  Assume  an  initial  flaw  size  for  the  primary  and  secondary  crack  lengths, 

Qp  and  a^. 

2.  Extend  the  primary  crack,  zlcTp  =  O.la^. 

3.  Calculate  the  stress-intensity  factors  (SIF)  for  both  the  primary  and 
secondary  cracks  and 

4.  Calculate  the  crack  growth  rate  of  the  primary  crack,  dOp/dN,  using 
Walker's  equation. 

5.  Calculate  the  cycle  interval  corresponding  to  the  primary  crack 
extension,  AN=  (dN/dap)Aap. 

6.  Extend  the  secondary  crack  over  the  cycle  interval,  Aa^  =  (da/dN)AN. 

7.  Sum  the  crack  lengths  for  the  primaiy  and  secondary  cracks,  a  =  a  + 

Aa. 

8.  Calculate  the  residual  strength  for  the  primary  and  secondary  cracks. 

Steps  2  through  8  are  repeated  until  the  second  breakthrough  occurs. 

REPRESENTATIVE  ORIGINAL  EQUIPMENT  MANUFACTURER  fROENf)  APPROACH. 

Results  from  the  damage  tolerance  analysis  methodology  in  RAPID  were  also  compared  to  results 
generated  using  a  Representative  Original  Equipment  Manufacturer  (ROEM)  approach.  In  the 
ROEM  approach,  the  method  used  to  determine  the  critical  fastener  location  and  the  fastener  load 
is  the  same  as  that  used  in  RAPID.  However,  the  initial  flaw  crack  geometry  and  the  crack 
growth  simulation  procedure  differ  from  those  in  RAPID  shown  in  figure  1.  It  is  assumed  in  the 
ROEM  approach  that  a  primary  and  a  secondary  cracks  of  sizes  0.05"  and  0.005",  respectively, 
exists  at  both  sides  of  the  critical  fastener  hole.  In  addition,  it  is  assumed  that  secondaiy  cracks  of 
size  0.005"  exists  on  both  sides  of  the  remaining  fastener  holes.  The  growth  of  each  crack  is 
assumed  to  occur  simultaneously  but  independently.  That  is,  the  interaction  between  the  crack 
tips  is  ignored  in  the  stress-intensity  factor  solutions  for  the  individual  cracks. 

In  the  ROEM  approach,  three  examples  representing  typical  fuselage  skin  repair  configurations 
were  analyzed.  Repair  type  I  consists  of  a  single  external  doubler  mechanically  fastened  over  a 
cutout  using  two  types  of  fasteners  as  shown  in  figure  7.  Repair  type  11  consists  of  an  internal  and 
external  doubler  mechanically  fastened  using  three  fastener  types  as  shown  in  figure  8.  Repair 
type  III  consists  of  a  skin  cutout  and  two  external  doublers  mechanically  fastened  using  three 
fastener  types  as  shown  in  figure  9.  The  crack  growth  and  inspection  schedule  were  calculated 
using  RAPID  and  the  ROEM  approach. 
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Doubler: 

2024-T3  aad  Sheet 
Thickness  =  0.07r 
Rivets: 

0  NAS1097-E6 
•  HL326-6 
Pitch  =i.(r 
Edge  distance  =  0.5* 

FIGURE  7.  DEFINITION  OF  REPAIR  TYPE  I  USED  IN  ROEM  ANALYSIS 
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Skin: 

2024-T3  aad  Sheet 
Thickness  =  0.063" 

Doubler: 

2024-T3  aad  Sheet 
Thickness  =  I 

Rivets:  (inner) 

0  NAS1097-E6 
•  HL326.6 
0)  MS20470-AD5 
Pitch  =1.0" 

Edge  Distance  =  0.5" 


FIGURE  9.  DEFINITION  OF  REPAIR  TYPE  IH  USED  IN  ROEM  ANALYSIS 


RESULTS  AND  DISCUSSION 


FASTENER  LOADS. 


A  comparison  of  the  calculated  fastener  loads  using  RAPID  and  FRANC2DL  is  shown  in  figure 
10.  As  shown  in  the  figure,  the  model  consists  of  a  skin  and  doubler  strips  fastened  using  three 
fasteners.  On  one  end,  the  doubler  is  held  fiixed,  and  on  the  other  end,  a  unit  stress  is  applied  to 
the  skin.  These  boundary  conditions  model  the  skin  cutout.  As  indicated  in  this  figure,  the 
critical  fastener  with  the  highest  load  calculated  using  FRANC2DL  was  approximately  1  percent 
lower  then  that  calculated  using  RAPID. 

In  another  comparison  of  the  fastener  load  calculations,  a  different  model  was  considered  as 
shown  in  figure  11.  As  shown  in  the  figure,  the  model  consists  of  a  skin  and  doubler  strips 
fastened  using  five  fasteners.  On  one  end,  both  the  skin  and  doubler  are  held  fixed,  and  on  the 
other  end,  a  unit  stress  is  applied  to  the  skin.  In  this  figure,  results  published  by  Swift  [1]  was 
used  as  a  baseline  solution  to  ensure  that  the  finite  element  model  used  in  the  FRANC2DL 
analysis  yielded  sensible  results.  As  indicated  in  the  figure,  the  fastener  loads  calculated  using 
FRANC2DL  were  in  good  agreement  with  results  published  by  Swift  [1].  For  the  critical  fastener 
with  the  highest  load,  the  difference  in  the  results  was  less  than  1  percent,  figure  11. 
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Fastener  Number 

FIGURE  10.  CALCULATED  FASTENER  LOADS  FROM  RAPID  AND  FRANC2DL 


Fastener  Number 

FIGURE  1 1 .  FASTENER  LOADS  CALCULATED  USING  FRANC2DL  AND  RESULTS 

FROM  REFERENCE  1 
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STRESS-INTENSITY  FACTOR  SOLUTIONS. 

Using  FRANC2DL,  the  SIF  as  a  function  of  crack  length  for  the  primary  crack  was  calculated  for 
all  three  cases  and  compared  with  RAPID  calculations  as  shown  in  figure  12.  In  general,  the  SIF 
solutions  calculated  for  Case  1  were  the  largest  values  while  the  solutions  calculated  using  RAPID 
were  the  lowest  values.  Results  for  Cases  2  and  3  were  in  the  middle.  The  trends  for  all  results 
were  the  same;  the  SIF  values  increased  with  crack  length.  As  the  crack  tip  approached  the  first 
hole,  the  SF  values  increased  substantially.  The  local  net  section  stresses  increase  as  the  crack  tip 
approaches  the  hole.  The  concentration  of  stress  along  the  hole  boundary  causes  the  SF  to 
increase.  A  step-wise  decrease  in  the  value  of  the  SF  occurred  immediately  after  the  first 
breakthrough.  Lower  stresses  in  the  net  section  ligament  length  between  the  crack  tip  and  the 
second  fastener  hole  caused  a  decrease  in  the  value  of  the  SF.  The  SF  values  then  increased 
with  increase  in  crack  length.  The  rate  of  increase  of  the  SF  decreased  as  the  crack  length 
increased.  As  the  crack  approached  the  second  hole,  the  rate  of  increase  of  the  SF  again 
increased  due  to  the  hole  boundary  effect. 
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FIGURE  12.  STRESS-INTENSITY  FACTOR  AS  A  FUNCTION  OF  CRACK  LENGTH 

For  crack  growth  within  two  hole  diameters  of  the  initial  crack,  there  was  good  agreement 
between  the  RAPID  results  and  all  three  FRANC2DL  results,  with  a  maximum  difference  of  2 
percent.  Modeling  the  loaded  hole  in  front  of  the  crack  tip  only  marginally  increased  the  SF 
solution.  As  the  crack  length  increased  and  approached  the  first  hole,  the  difference  in  values 
increased  substantially.  Here  the  effect  of  modeling  a  loaded  hole  in  front  of  the  crack  tip  can  be 
seen  by  comparing  results  from  Case  1  with  the  other  results.  The  SF  results  for  Case  1  were 


o  RAPID  Analysis 
□  FRANC2DL  Analysis  (Case  #1) 
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substantially  higher.  In  addition,  modeling  an  array  of  holes  in  front  of  the  crack  tip  did  not  affect 
the  SIF  results  as  seen  by  comparing  results  of  Cases  2  and  3.  Thus,  modeling  an  array  of  holes  in 
front  of  the  crack  tip  is  not  necessary.  The  RAPID  approximation  does  not  fiilly  capture  the 
singularity  as  the  crack  tip  approaches  the  hole,  however,  this  occurs  over  a  very  short  crack 
extension  length. 

After  the  first  breakthrough  condition,  there  was  a  maximum  of  a  5  percent  difference  in  the 
solutions  for  crack  lengths  up  to  approximately  25  percent  of  the  ligament  length  between  the  first 
and  second  breakthrough  holes.  From  this  point  forward,  the  difference  increased  as  the  crack 
length  increased.  Contributing  to  this  increased  difference  is  the  effect  due  to  the  finite  width  of 
the  finite  element  models  (RAPID  assumes  crack  growth  in  an  infinite  body).  For  crack  lengths 
over  50  percent  of  the  ligament  length  between  the  first  and  second  breakthrough  holes,  the  tip- 
to-tip  a/iv  ratio  was  approximately  0.2.  For  this  a/W  ratio,  it  is  expected  that  the  SIF  calculated 
using  FRANC2DL  would  be  5  percent  higher  than  the  RAPID  results  due  to  the  finite  width 
effect. 

CRACK  GROWTH. 


The  length  of  the  primary  crack  as  a  function  of  the  number  of  flights  was  calculated  using 
FRANC2DL  and  compared  with  RAPID  results  as  shown  in  figure  13.  In  this  figure,  the  primary 
crack  length,  a^,  is  measured  from  the  tip  of  the  crack  to  the  centerline  of  the  initial  hole.  In 
general,  the  results  are  in  good  agreement.  Differences  in  the  crack  growth  predictions  are  due 
mainly  to  the  differences  in  the  SIF  calculations  as  discussed  previously.  For  crack  lengths  up  to 
approximately  50  percent  of  the  first  breakthrough  condition,  the  corresponding  number  of  flights 


FIGURE  13.  CRACK  LENGTH  AS  A  FUNCTION  OF  THE  NUMBER  OF  FLIGHTS 
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calculated  using  RAPID  was  5  percent  higher  than  the  FRANC2DL  results.  As  the  crack  length 
increased  and  approached  the  first  hole,  the  number  of  flights  calculated  using  RAPID  was 
approximately  8  percent  higher  than  the  FRANC2DL  result.  Up  to  the  second  breakthrough,  the 
number  of  flights  calculated  using  RAPID  was  10  percent  higher  than  the  FRANC2DL  result. 
Better  agreement  would  be  obtained  using  Cases  2  and  3  which  better  resemble  the  loading 
conditions  used  in  RAPE). 

The  crack  length  as  a  function  of  the  number  of  flights  was  calculated  using  the  ROEM  approach 
for  the  three  repair  types  and  compared  with  results  generated  using  RAPID  as  shown  in  figures 
14-16.  In  general,  a  slow  crack  growth  process  was  predicted  up  to  the  first  breakthrough 
condition  for  both  analysis.  At  the  first  breakthrough  condition,  the  ROEM  analysis  predicts 
sudden  catastrophic  fi'acture  of  the  repair  skin.  The  RAPED  analysis  predicts  a  more  stable  crack 
growth  process  where  at  least  two  breakthrough  conditions  occurred  prior  to  catastrophic 
fi’acture  of  the  rep^  skin.  The  difference  in  the  results  using  the  two  methods  is  due  to  the 
differences  in  the  initial  flaw  assumptions.  In  the  ROEM  approach,  primary  and  secondary  cracks 
of  lengths  0.05"  and  0.005",  respectively,  exists  at  both  sides  of  the  critical  fastener  hole  and 
secondary  cracks  of  length  0.005"  exists  on  both  sides  of  the  remaining  fastener  holes.  The 
growth  of  each  crack  is  assumed  to  occur  simultaneously  but  independently.  In  the  RAPED 
analysis,  the  Case  2  crack  growth  scenario  shown  in  figure  1  was  assumed. 

The  crack  growth  results  for  repair  type  I  are  shown  in  figure  14.  For  the  ROEM  analysis,  the 
first  breakthrough  condition  occurred  at  47,000  flights,  where  the  initial  primaiy  crack  grew  to 
the 


N  (flights) 

FIGURE  14.  TOTAL  CRACK  LENGTH  AS  A  FUNCTION  OF  THE  NUMBER  OF  FLIGHTS 

FOR  TYPE  I  REPAIR 
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adjacent  hole  and  the  initial  secondary  cracks  grew  independently  to  lengths  of  0.525".  At  these 
crack  growth  lengths,  all  cracks  linked  up  to  form  a  single  crack  resulting  in  sudden  catastrophic 
fracture  of  the  repair  skin.  For  the  RAPID  analysis  shown  in  figure  14,  the  first  breakthrough 
occurred  at  41,000  flights  where  the  initial  primary  crack  grew  to  the  adjacent  hole.  Secondary 
cracks  then  initiated  at  opposite  sides  of  the  holes  as  shown  schematically  for  Case  2  in  figure  1. 
The  subsequent  crack  growth  rate  increased  with  the  second  breakthrough  occurring  at  45,515 
flights.  The  RAPID  analysis  indicated  that  catastrophic  fracture  of  the  repair  skin  occurred  after 
the  third  breakthrough  at  46,805  flights. 

Trends  for  the  results  for  repair  types  n  and  HI  shown  in  figures  15  and  16,  respectively,  are 
similar  to  those  in  figure  14.  For  the  ROEM  analysis,  the  first  breakthrough  condition  for  repair 
type  n  and  HI  occurred  at  53,100  and  53,586  flights,  as  shown  in  figures  15  and  16.  After  the 
first  breakthrough  condition,  sudden  catastrophic  fracture  of  the  repaired  skin  was  predicted  using 
the  ROEM  analysis.  For  the  RAPID  analysis,  the  first  breakthrough  condition  for  repair  type  11 
and  in  occurred  at  48,416  and  48,928  flights,  as  shown  in  figures  15  and  16.  Fracture  of  the 
repaired  skin  after  the  second  breakthrough  condition  was  predicted  using  the  RAPED  analysis  at 
53,100  and  53,586  flights  for  repair  types  H  and  m,  respectively. 
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HGURE  15.  TOTAL  CRACK  LENGTH  AS  A  FUNCTION  OF  THE  NUMBER  OF  FLIGHTS 

FOR  TYPE  n  REPAIR 
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N  (flights) 


FIGURE  16.  TOTAL  CRACK  LENGTH  AS  A  FUNCTION  OF  THE  NUMBER  OF  FLIGHTS 

FOR  TYPE  III  REPAIR 


RESIDUAL  STRENGTH. 

Residual  strength  as  a  function  of  crack  length  for  the  primary  crack  was  calculated  using  RAPID 
and  FRANC2DL  as  shown  in  figure  17.  The  residual  strength,  S^s,  was  calculated  using  equation 
5.  Case  1  shown  in  figure  5  was  used  in  the  FRANC2DL  analysis.  In  general,  the  trends  for  the 
results  were  the  same;  the  residual  strength  values  decreased  with  crack  length.  As  the  crack  tip 
approached  the  first  hole,  the  residual  strength  values  decreased  substantially  due  to  hole 
boundary  effects.  The  residual  strength  is  inversely  proportional  to  the  SIF  term  (denominator)  in 
equation  5.  As  the  crack  tip  approaches  the  hole  boundary,  the  local  net  section  stresses  increase 
causing  the  SIF  to  increase  and,  thus,  the  residual  strength  to  decrease.  A  step-wise  increase  in 
the  value  of  the  residual  strength  occurred  inunediately  after  the  first  breakthrough  due  to  the 
decrease  in  SIF.  The  stresses  in  net  section  ligament  between  the  crack  tip  and  the  second 
fastener  hole  deceased  resulting  in  the  decreased  value  of  the  SIF.  With  further  increase  in  crack 
length,  the  residual  strength  decreased.  As  the  crack  approached  the  second  hole  the  residual 
strength  again  decreased  due  to  the  hole  boundary  effect. 
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FIGURE  17.  RESIDUAL  STRENGTH  AS  A  FUNCTION  OF  CRACK  LENGTH 


The  differences  in  the  residual  strength  shown  in  figure  17  are  due  mainly  to  the  differences  in  the 
SIF  calculated  using  RAPID  and  FRANC2DL.  The  largest  differences  occurred  as  the  crack 
length  approached  the  first  and  second  breakthrough  conditions. 


INSPECTION  INTERVAL. 


The  inspection  interval  as  a  function  of  the  detectable  crack  length  was  determined  using  the 
ROEM  approach  for  three  repair  types  and  compared  with  results  generated  using  RAPID  as 
shown  in  figures  18-20.  The  inspection  interval,  was  calculated  using  equation  9.  The  first 
data  point  is  the  inspection  threshold.  Nth,  defined  using  equation  6.  As  shown  in  figures  18  -  20, 
excellent  agreement  was  obtained  between  the  RAPID  and  ROEM  results  up  to  the  first 
breakthrough  condition.  The  discrepancies  in  the  results  using  the  two  methods  is  due  to  the 
differences  in  the  initial  flaw  assumptions.  In  the  ROEM  approach,  primary  and  secondary  cracks 
of  lengths  0.05"  and  0.005",  respectively,  exists  at  both  sides  of  the  critical  fastener  hole  and 
secondary  cracks  of  length  0.005"  exists  on  both  sides  of  the  remaining  fastener  holes.  The 
growth  of  each  crack  is  assumed  to  occur  simultaneously  but  independently.  In  the  RAPID 
analysis,  the  Case  2  crack  growth  scenario  shown  in  figure  1  was  assumed. 
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FIGURE  20.  INSPECTION  INTERVAL  FOR  REPAIR  TYPE  HI 


CONCLUDING  REMARKS 

In  this  study,  the  simplified  damage  tolerance  analysis  methodology  in  the  newly  developed  rep^ 
software.  Repair  Assessment  Procedure  and  Integrated  Design  (RAPID),  was  presented  and 
evaluated.  Results  generated  using  RAPID  compared  well  Avith  results  generated  using  a 
Representative  Original  Equipment  Manufacturer  (ROEM)  method  and  a  special  purpose  finite 
element  program,  the  Fracture  Analysis  Code  for  2-Dimensional  Layered  (FRANC2DL) 
structures.  Comparisons  were  made  in  terms  of  the  fastener  loads,  stress-intensity  factor  (SIF) 
solutions,  crack  growth,  residual  strength,  and  inspection  intervals. 

The  fastener  loads  calculated  using  RAPID  and  FRANC2DL  were  in  good  agreement.  There  was 
a  1  percent  difference  in  the  value  of  the  load  at  the  critical  fastener  hole. 

Trends  of  the  SBF  results  calculated  using  RAPID  and  FRANC2DL  were  the  same;  the  SIF  values 
increased  vdth  crack  length.  As  the  crack  tip  approached  the  first  hole,  the  SIF  values  increased 
substantially.  The  local  net  section  stresses  increase  as  the  crack  tip  approaches  the  hole.  The 
concentration  of  stress  along  the  hole  boundary  causes  the  SIF  to  increase.  A  step-wise  decrease 
in  the  value  of  the  SIF  occurred  immediately  after  the  first  breakthrough.  Lower  stresses  in  the 
net  section  ligament  length  between  the  crack  tip  and  the  second  fastener  hole  caused  a  decrease 
in  the  value  of  the  SIF.  The  SBF  values  then  increased  with  increase  in  crack  length.  There  were 
noticeable  differences  in  the  results  as  the  crack  tip  approached  a  hole  boundary;  the  values  of  the 
SIF  calculated  using  FRANC2DL  were  higher  than  those  calculated  using  RAPID.  However,  this 
occurred  over  a  short  crack  extension  length. 
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The  crack  length  as  a  function  of  the  number  of  flights  calculated  using  RAPID  was  compared 
with  results  using  FRANC2DL.  Good  agreement  was  obtained  between  the  RAPID  and 
FRANC2DL  results  with  a  8  percent  difference  in  the  number  of  flights  to  reach  the  first 
breakthrough  condition.  The  difference  in  the  results  was  due  to  differences  in  the  SIF, 

Both  RAPID  and  the  ROEM  approach  were  used  to  calculate  the  crack  length  as  a  function  of  the 
number  of  flights  for  three  repair  types.  In  general,  a  slow  crack  growth  process  was  predicted  up 
to  the  first  breakthrough  condition  for  both  analysis.  At  the  first  breakthrough  condition,  the 
ROEM  analysis  predicts  sudden  catastrophic  fracture  of  the  repair  skin.  The  RAPID  analysis 
predicts  a  more  stable  crack  growth  process  where  at  least  two  breakthrough  conditions  occurred 
prior  to  catastrophic  fi’acture  of  the  repair  skin.  The  difference  in  the  results  using  the  two 
methods  is  due  to  the  differences  in  the  initial  flaw  assumptions. 

The  residual  strength  calculations  made  using  RAPID  and  FRANC2DL  were  in  good  agreement, 
with  the  largest  differences  occurring  as  the  crack  tip  approached  the  fastener  hole.  The 
difference  in  the  results  was  due  to  differences  in  the  SIF. 

Finally,  the  inspection  interval  calculated  using  RAPID  for  the  three  different  repair  types  were  in 
good  agreement  with  results  generated  using  the  ROEM  approach  up  to  the  first  breakthrough 
condition.  The  difference  in  the  results  using  the  two  methods  is  due  to  the  differences  in  the 
initial  flaw  assumptions. 
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